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Abstract. An asymmetric exclusion process comprising positive particles, negative 
particles and vacancies is introduced. The model is defined on a ring and the dynamics 
does not conserve the number of particles. We solve the steady state exactly and show 
that it can exhibit a continuous phase transition in which the density of vacancies 
decreases to zero. The model has no absorbing state and furnishes an example of a 
one-dimensional phase transition in a homogeneous non-conserving system which does 
not belong to the absorbing state universality classes. 
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One dimensional driven diffusive systems have been a subject of extensive study 
in recent years [JXJ, 0, [|. A prototypical model for studying these systems is the 
asymmetric exclusion process in which particles move stochastically with a preferred 
direction and hard-core exclusion J3|. In this model the local dynamics conserves the 
particles. Asymmetric exclusion processes with a single species of particles do not 
exhibit phase transitions on a ring geometry. On the other hand open systems in which 
the particle number is not conserved at the boundaries can exhibit a variety of phase 
transitions f||, [5|]. Generalisations of these models to more than one species of particles 
have shown that phase transitions and long-range order do exist in these systems even 
on a ring geometry f|, [7], || |9|, . 

Generally speaking one-dimensional models that exhibit phase transitions are either 
(a) of the asymmetric-exclusion-process type with a drive and conserving bulk dynamics 
or (b) have non-conserving bulk dynamics with one or more absorbing states. The 
latter case is related to the directed percolation or other absorbing state universality 
classes [[II]. It is therefore of interest to find a one-dimensional model where the bulk 
dynamics is not conserving that exhibits a phase transition unrelated to the absorbing 
state universality classes. 

In this letter we introduce such a model and solve for its steady state exactly. The 
model is defined on a ring of L sites where each site can be occupied by either a positive 
(+) particle, a negative (— ) particle or a vacancy (0). The model evolves through the 
following conserving rates: 

+ 0^0+ 0-^-0 +-^-+ (1) 

i 

augmented by the following non-conserving rates 

+ 0^ 00 -^00. (2) 

li v ' 



Thus, the model generalizes the model of Derrida et al [?J and Arndt et al [|1(J through 
the introduction of the process of creation and annihilation of particles. Note that any 
configuration can be reached from any other except states with no vacancies which are 
not dynamically accessible. 

The rate w controls the density of vacancies, denoted by 9, in the system. For large 
w one expects the density of vacancies to be high and for small w it is expected to be 
low. We will show that for q < 1 the dependence of 9 on w is not smooth and that for 
small enough w , 9 is zero in the thermodynamic limit. Thus the model exhibits a phase 
transition from a "fluid" phase with a finite 9 to a "maximal current" phase where 9 
is zero (the nomenclature will be explained below). The transition occurs at a critical 
value w c > 0. The phase transition is found to be continuous with 

0~\w-w c f, (3) 

where (3=1. For q > 1, on the other hand, the system is always strongly phase 
separated f|, El 03 with a single zero and two extensive pure domains of positive and 
negative particles. Here there is no phase transition as w is varied. 



Phase transition in a non-conserving driven diffusive system 



3 



A powerful technique for solving the steady-states of asymmetric exclusion processes 
is the matrix product ansatz [I^j. This involves representing the steady-state weights as 
the trace of a product of matrices which depends on the microscopic configuration. The 
matrices then obey certain algebraic rules which are derived from the dynamics of the 
model. This technique has yielded the exact behaviour at the various phase transitions 
in many asymmetric exclusion models. The transitions solved within the matrix product 



have been found to be robust for a large class of systems |13fl . However, in most models 



solved so far using the matrix product, particle numbers are conserved in the bulk (two 
exceptions are fl4| and fl5|]). In this work we employ the matrix product technique to 
obtain an exact solution for the steady state and phase transitions of the model defined 
by It turns out that despite the non-conserving dynamics of the present model 

we can use the same matrices that have been previously used to solve models with 
conserving dynamics. 

We now proceed to outline the matrix product solution for the steady state. The 
steady state weight of a configuration C is represented by the trace of a product of 
matrices: 

L 

W L {C)=TrH[6 Ti!+ D + 6 Ti> -E + 6 nj0 A] , (4) 
i=i 

where L is the system size and Tj = +, — , if site j is occupied by a + , — or respectively. 
That is, a matrix D (E) represents a positive (negative) particle and A represents a 
vacancy. It is easy to show using the technique of |12[ [?]] that D, E and A should satisfy 

DE - qED = D + E ■ DA = AE = A ; AA = wA (5) 

to give the correct steady-state weights. These equations can be satisfied by choosing 

A = w\V)(V\, (6) 

where (V\V) = 1, D\V) = \V) and (V\E = (V\. Then D, E and (V\ are identical to the 
matrices and vectors studied in Jl6|, [L7j where they are used to solve the single-species 
partially asymmetric exclusion process with open boundaries. 

We now wish to calculate the normalisation, i.e. the partition function, which is 
given by the sum of the weights (J|) of all accessible configurations: 

Z L = Tr [(A + D + E) L - (D + E) L ] (7) 

Using the form (|6|) of A we first write the sum of the weights of configurations with 
exactly M vacancies on a lattice of size L: 

M oo 

Zlm = » m I1E (V\C^\V) S e ^ l _ m (8) 
m=i »v=o 

where C = D + E. The sum over each n M corresponds to the possible number of particles 
between two consecutive zeros. The delta function enforces the constraint that the total 
number of particles equals L — M and the factor w M arises from the M zeros. This 
form neglects the degeneracy in placing a given configuration {n^} on a ring geometry, 
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which is bounded from above by L. It is straightforward to check using bounds on the 
true partition function that this does not affect any of the results present here. 

In order to calculate the partition function Zl = Em=i %l,m it is convenient to 
replace the delta function by a contour integral. Taking the limit of the sum over M to 
infinity yields 

00 r dz w M M 00 
Zl= E /tE^mTT n E z^(V\C^\V) (9) 

(10) 



7 L-M+1 



dz 



2m z L+1 
where U(z) is defined as 



zwU (z) 
— zwU(z) 



U(z) = Y J z n G n with G n = (V\C n \V). (11) 

n=0 

The weight G n has been studied before. It is the normalisation sum of the single species 
partially asymmetric exclusion model on an open lattice of size n and with particle 
injection rate 1 at the left boundary and removal rate 1 at the right boundary [[U| |T^| . 
Below, we shall consider the density of vacancies given by 

M ,. wd\nZ L . . 

6= lim — = lim — — - , (12) 

where M is the average number of vacancies in the system. We shall also consider 
the current Jl of positive particles (which is equal to that of the negative particles). 
We define the current in the same way as for the conserving system although here the 
density is not conserved. Thus taking into account the inaccessibility of configurations 
with no vacancies we find 

J L = — Tr \{DA + DE- qED) (C + A) L ~ 2 - (DE - qED) C L ~ 2 ' 

= %± • (13) 
Zl 

where we have used the algebraic rules (||) and the form (0) of Z L . 
We now discuss the two distinct cases q < 1 and q > 1 separately: 



(i) The case of q < 1 

Here the normalisation Zl can be evaluated for large L from the integral (10) by 
the saddle point method. This method amounts to working in the grand canonical 
ensemble. The term in the square brackets of (|T0"| ) is then just the grand canonical 
partition function and Zl ~ (z*) ~ L where z* is the saddle point value of z, i.e. the 
fugacity. Thus, in the thermodynamic limit, we see from (p~3| ) that Jl — > z* = J which 
identifies the particle current in the system as the fugacity. Also, we see from (|T2| ) that 
the density of vacancies is given by 
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Thus one may identify 9 as the order parameter and — In z* = 
the free energy density. 

The saddle point equation for the integral ( [HDD may be written as 



In J as the analogue of 



L 



wzU (z) 



zU'(z) 
U[z) 



+ wzU(z) 



(15) 



For L — > oo this equation is satisfied either by 1 — wzU(z) ~ 0(1/ L) or by the term 
inside the brackets being of order L. To analyse which of the two scenarios pertains 
one needs to consider the properties of the increasing function U(z). Note from (pT|) 
that the convergence of U(z) is determined by the large n form of G n . For q < 1 this 

[L7|] to behave for large n as 

K -n 



quantity is known |16 
G ~ 



n 



3/2 



where 



K 



1-q 



and a = —= 



n 



i\3 



(16) 



(17) 



From fll6|) one deduces that ?7(z) converges for z < K and diverges as z — > K. 

Thus, when 1 — wzll(z) = for some z < K we solve the saddle-point equation ([TJD by 
choosing z* so that 1 — wz*U(z*) ~ 0(1/ L). On the other hand when 1 — wzU(z) > 
for all z < K, the divergence must come from the term in the square brackets of 
equation flT5| ) and we need = K(l — C(l/L 2 )). This can be deduced by noting that 
U'(z) ~ \K — z\~ l l 2 for z — > X. Thus, in the thermodynamic limit L 
to as w is decreased to iy c , given by 

1 



oo, 2; increases 



KU(K) ' v ' 

Any further decrease in w leaves the value of z* unchanged. Using equation flTjj ) one 
can see that for w > w c , 9 decreases as w decreases while for w < w c , 9 = 0. Using 



U(z) - U(K) 



\K 



1 1/2 



it is easy to show by expanding ([15) that 9 ~ \w — w c \ as 



w 



w c , recovering 



Since z* is the current of particles and it saturates at K = (1 — q)/A for w < w c we 
refer to this phase as the maximal current phase. We refer to the phase w > w c as the 
fluid since the typical configuration is a disordered arrangement of +, — and Os. 



(ii) The case of q > 1 

For q > 1 it is known [17] that the n dependence of G n is given by 



G n 



Q 



q 



n 2 /4 



(19) 



In this case U(z) diverges for all z > and we have to impose cutoffs on the sums in 
equation @. That is, we write 



Z, 



E 

M=l 



dz w 



M 



fh-M 



M 



2ni z 



L-M+l 



(20) 



. n=0 
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From ( |T9D it is clear that the dominant contribution to the sum is when M = 1 and 
n — L — 1. Then Ji = Z^^/Z^ ~ (q — 1) q~ L l 2 ~ l . Thus, the current of particles 
is exponentially small in the system size and the density of vacancies is zero. This 
corresponds to a strongly phase separated state with a single vacancy followed by an 
extensive block of positive particles (to the right of the vacancy) followed by an exten- 
sive block of negative particles (to the left of the vacancy). 



Finally, we discuss the transition between q < 1 and q > 1. Using equations ([11]) 
and ( |TTD it is easy to show that as q — > 1 from below, KU(K) tends to zero. Thus using 
equation (|T^) we see that w c diverges as q — > 1. Therefore at the q = 1 transition the 
system changes from a maximal current phase J = (1 — q)/A for q < 1 to a strongly 
phase separated state with J = 0(q~ L / 2 ) for q > 1. In both of these phases the density 
of vacancies is zero. Therefore the transition has no relation to the non-conserving 



dynamics and is instead related to the reversal of the bias |17 



It is interesting to make a comparison between the transition from the fluid to the 
maximal current phase and the denaturation transition in DNA where the two strands 
of the molecule unbind at a certain temperature. In models of this transition [I3L Il9| one 



assigns a Boltzmann weight w to bound base pairs and a weight K~ n /n c to unbound 
segments of length n. Here K is a non-universal constant whereas c is a universal 
exponent depending on the dimension and self-avoidance properties of the unbound 



DNA. For example, using a random walk model yields c = 3/2 in three dimensions [18 



As temperature is raised and w decreases there is an unbinding transition where the 
fraction of bound base pairs, 8, vanishes. The DNA models and the present model 
can be related by identifying vacancies with bound pairs and blocks of particles with 
unbound DNA loops. Then G n corresponds to the weight of DNA loop of length n and 
the grand canonical partition functions of the two systems are the same . 

The difference between the systems is that in the DNA the exponent c = 3/2 is a 
result of the fact that the one-dimensional molecule is embedded in three dimensions. 
However, for the model considered here the exponent 3/2 arises from truly one- 
dimensional phenomena related to the currents flowing in the system. This can be seen 
from the relation between G n and the current of the single-species partially asymmetric 
exclusion process. More explicitly, for such a process the current in a system of size n is 
given by G n -i/G n . Thus, the n dependence of G n is related to the n dependence of the 
current in a single species system of size n. In the context of the two species model we 
have studied here such a system corresponds to an uninterrupted block of particles, of 
length n, bounded between two vacancies. This picture has recently been used to study 
the conditions for phase separation in conserving models |2T[ . 

The model can easily be generalised to contain two additional parameters a and (3 
by modifying the rates ([[]) to read 

+ 0^0+ O-A-0 +-^-+ (21) 
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and the rates (|2]) to 

3w aw 

+ 0^00 0-^00. 22 

11 v ' 

Then one can take A = w\V)(W\ where (3D\V) = \V) and a(W\E = (W\. This 
generalisation should allow for a richer phase diagram than that presented here. 

Finally, we point out that the transition from the fluid to the maximal current 
phase is lost if particles are not conserved inside a block consisting only of particles. 

For example, a process where H > regardless of its neighbours destroys the maximal 

current phase. 
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